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Thermal Analysis of Antenna Backup Structure 
Part I Methoclology Development 

D. Schonfetd and F. L. Lansing 
OSN Engintvrino S«ctton 


An analvtic method is devised for predicting the temperature disniiK.non in typkni 
onteniKt s:ructurai back-up members. The remits are in good agreement with those 
obtained by a numerical shifoting method The analytic method developed in this work 
has shown a good potential in greatly simplifying the thermal analysis process for 
complex back-up antenna structures 


I. Introduction 

Future requirements for dec^ space navigation and telem- 
etry link communications point to the need for higher 
frequenc) bands. At these higlier frequencies, data rates are 
better but antenna surface accuracy requirements are more 
severe than those of the presently used S- or X-bands. Guided 
bv the ac».urac\ needs, studies are lum in pn>gre« on the 
feasibilit) of constructing large Ka-ban I antennas which will 
use frequencies at 32 GH/ or higher (Ref. 1 ). 

T*ie construction requirements for these Ka-band antennas 
are much more demanding than those for presently used 
devices. For example, the reflective surfaces must be pv>si- 
(loned 4iid aligned much more accurately and must be 
maintained during the life of the antenna against varying 
gravity, wind, and thermal loadings. Fnvironmenlal changes 
such as temperature variations can have deleterious effects on 
the antenna performance. Field thermal measuremeh;.> 
(Ret. 2> and analytic investigations are being conducted to 
study the environmental effects on the structural members of 
the large b4-m antenna. In the analytical investigation, which 
IS the subiect of this irticle. we are interested in simulating the 
temperature distribut:i>n throughout the complex back-up 


structural members given their phy^sical properties, geomethc 
arrangement, and environmental conditions, such as air tem- 
perature. wind veliv'ity. and solar inadiation (insolation). 

The simulation of the temperature pattern of a structural 
member can be done by ^conventional” numerical finite 
difference methods in heat transfer which divide the member 
under a>nsideration into nodes and then apply heat oalance 
equations to each node. Rather than solving for an excessively 
large network of nodes F a complex antenna back-up 
structure, a new method is developed in this paper to save 
effort and time. The method is analytical in nature and relies 
i>n deriving a universal relationship for the temperature 
varutions and heat fluxes within each member. The methodol- 
ogy is described only in this article and will be followed by 
additional applications iii subsequent TDA reports. 

It. Methodology Development 

Consider a simple bar, of length /.,as sketdied in Fig. I 
subjected to solar radiation, cmduction. convection, and 
radiation heat exchange with the ambient atr. The cross 
section of the bar. .4. and its material properties are assumed 
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CQostaat akng the bar’s kagth. For an ekment of iDaemental where F, stands for the tea ipe ratme, 7, whSe 7] and 7] for 

length dx akmg the bar axis, whose temperature is T{X), the dTIdx and d^Tfdic^ re s pe c thw l y. a^ the supeisci^ (') for 


oonvectioo losses am gwen by: 


d/dx-lbe constants C,,C 3 andCj aredetfaiedas: 


dq^ = kPdx(T-T^ 

0) 

c 

« kA 

(10) 

the ndietne losses by: 




dq^ = eo^Pdx{T*-Tp 

(2) 

C 

* kA 

(II) 

the absorbed pctoon of the solar irradiation by: 


- _ kP _ €OP^ _4 oZP 
~ ~kA V-tt* 

(12) 


dq^ * alPdx (3) 


This numerical method wffl serve as a useful check on the 
and the Jiange in conductive heat transfer stored in the analytical results to be devdoped fatter, 
element b>: 



-) 
dx I 


dx 


(4) 


The complete heat balance equation at $teady>state conditions 
can be written as: 


d^T 

dx* 


• — (7- T)* (T* - 

kA ^ V kA kA 


(5) 


This analytic method of stdvmg Eq. (S) using perturbatiao 
requirm that the equation must become linear in the limit of 
some selected small par a meter. The pananeter C, in Eq. (10) 
is small because of the Stefan^oltzma constat o of the 
order 10^'*, while C] and from Eqs. (11) and (12) are of 
tlie order 1 and greater. The perturbation method assumes that 
the solution can be written in the form: 


TXx) = 7^(x) + C, 7,(x) + C^7j(x) + • • (13) 


Equation (S). together with the boundary conditions 
at X = 0 , r * 7(0) 

at X = I , 7 = 7(i) 


( 6 ) 

(7) 


Then, by inserting (13) into (S), one (Stains a set of linear 
equations such as: 



C,7„-C3=0 


(14) 


fonns a nonlinear boundary value problem for T\X). in 
general, such problems cannot be solved analytically and 
numerical methods must be employed (Refs. 3 and 4). Several 
of these methods are presented in the paragraphs whidi fdlow. 

A. 


dFr 

--CJ, 

dx* * ‘ * 0 


^-C,73-C3 = 7t 


An outline of the ''shooting" method and its use is given in 
Ref. 4. Using the subroutine described in that reference, 
Eq. (5) can be written as a system of equations in the form: 


which must be solved sequentially. 


(15) 

(16) 


y' = V 


Y* “ C 




The disadvantage of applying this method for our case is 
(8) that once an expression for is found by solving Eq. (14), 
subsequent equations will become very complicated due to the 
presence of terms such as 7^, 7^, etc. For this reason, the 
(^^ perturbation method was not examined further. 
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C. U nsar i z H o n llBi hod 

Equation (5) will be reduced to a closed form solution for 
by Unearizifig the radiation term in Eq. (2) by writing it 
as: 


The ratio from Eqs. (20) and (21) is: 


JL 

* 




oJ 


(23) 


= h^PdxiT-TJ (17) 

By comparing Eq. ( 1 7) with Eq. (2) one obtains: 

= €0 + rjxrt r ) (i8) 


The initial value of can be obtained by making the 
approximation T ^ Hiis assumption is valid in cases where 
one would not expect too much difference between the 
temperature of the structure and that of the ambient air. Note, 
however, that such an assumption cannot be used for a 
structure in space w^.^^e the temperature of the structure is 
much different from the surrounding space temperature. (7^ - 
0 K). As the temperature pattern is known, a modified value 
of k^ can be obtained at the average Hnk temperature. 

By using Eq. (17) in Eq. (5), one obtains the linearized 
differential equation: 


ill 

dX^ 


L kA I kA 


T ^ 


g/p l 
kA J 


(1^) 


Hence, the physical meaning of represents a balance 
between the incoming insolation and the rate of radiative md 
convective losses, all with respect to the ambient temperature, 
r. Therefore. is an ^equilibrium** temperature which wffl 
be denoted by In terms of Eq. (22) becomes 

(7(A) - rj - (7(0) - g cosh v'5i ^ 

r(jr) = sinh y/^x 

sinh y/^L 

♦[no)-rjcosh + (24) 


m. Nunwrical Exmiple 

To check the validit>' of Eq. (24) versus the more accurate 
shooting method, let us consider a numerical example. In this 
example we assume a steel bar having the following geometry: 

Perimeter. F = 0305m (1 ft) 

Length. iL = 0.91m (3 ft) 

Cross sectional area. .4 = 5.81 X 10'^ m^ (0.0625 ft^) 


and by making use of the abbreviations 


The two end temperatures are kept at: 


^ = 


1^1 )F 

r c 

kA 


( 20 ) 


r(0) = 3II K(560*R) 
7(A) = 322K(580'*R) 




(A *h)P 


kA 


T + 


aJF 

kA 


{ 2 \) 


while the ambient temperature is: 

= 294K(530®R) 


one obtains the solution of Eq. ( 19) as: 


T{x) 


j^7(A) - 4- j - [r(0) - |-J cosh >/^A 


sinh 


sixth y/^X 


^nO) - “I" j cosh \^x + (22) 


For these temperatures, and assuming typical Goldstone 
values for wind speed [16.1 mfhx (10 mph)] and insolation 
(800 W/m^), we will use the following heat transfer coeffi- 
cients: 

conductivity, k = 45.0 W/mK (26.0 BTU/hr ft®R) 

radiation heat transfer coefficient h =3.1 W/m^ K (0.55 
BTU/hr ft^R) 

convective heat transfer coefficient, h = 22.7 W/m^ K (4.0 
BTUlirft^^R) 
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emissivity, e = O.S 
view factor, .^=1.0 
absorptivity, a - 0.4 

insolation, / = 800 W/m^ (254 BTU/hr ft^) 

For these values, ^ = 30.1 (2.8 and T= 308 K 

(5S2^R). The results of the footing method and the 
linearization method in Eq. (24) are shown in Tables 1 and 2. 
A [dot of the results of Table 2 is made in Fig. 2. llie excellent 
agreement between the analytical and numerical results vali- 
dates the use of Eq. (24). Note that the temperature 
distribution within the bar need not be monotonic; in this 
particular ca.<L% the temperature within the bar reaches a point 
where it is low. r than either of the end temperatures. 


IV. The MultMink Problem 

Consider now a system of 3 bars linked together as 
illustrated in Fig. 3. The physical and geometric properties of 
each bni can be di^erent and so can their exposure to solar 
radiation. The temperatures at the junctions are assumed 
to be unknown. The problem is to find a simple way of 
determining these temperatures and then, by using Eq. (24). to 
obtain the temperature distribution throughout the links. 

The solution relies on the fact that, at each junction of two 
or more members, the heat flux balance at steady state is 
written as: 


In Eqs. (25) through (27), the square brackets ace subscripted 
by the Unk number, while tlx derivitive dT/dx is to be taken 
at the ^beginning** (0) or the "*end** (L) of the respective link. 
The assignment of (0) and (L) to each link is completely 
arbitrary. Note that, based on Eq. (24): 

(28) 

(^\ ^ _v g^ _ . j. ) _ ^ 

y^/o Sinh>/Ji 

(29) 

By substituting Eq. (28) and (29) into Eqs. (25) through (27), 
(Mie gets three equations and ^ unknowns, the utdcnowns 
being the temperatures 7X£)^. and TtO)^. (i = 1, 2, 3). However, 
three of these temperatures are redundant because of the 
junction conditions: 


r(0), » 7Ti)j 

(3C) 

r(0)3 = no, 

(31) 

no, “ no)j 

(32) 




I./.*. 


..=0 


Since the only heat transfer at the junction points is through 
conduction, for the case we are considering, the following 
equations are valid: 



We choose to solve for T(0)^. (/= 1, 2, 3) and to use the 
following abbreviations: 


3 . = kA,- 

f t t 


V^- 


sin 




(33) 


= e^r,.(coshNA5‘i,-i) (34) 


6 ^ = cosh (35) 


Then, the system Eqs. (25) through (27), by using Eqs. (30) 
througli (32) and (24), can be written in matrix form as. 
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The solutions of Eq. (36) are given by: 


7-(0), 


r(0). 


■((», + (jj) 


-0| ' 



-«3 


-«3 

6,+«3_ 




-6,. 6^ 




(0, fi,) 

-»3 


(0, +0j) 

6 ,+ 63 _ 



(37) 


(38) 


r(0)3 


■6, +63 


(0,+0,)' 


63 + 63 

(0,*03) 


-«3 

(0.^03), 



( 39 ) 


whci?: 


6. +6, 



I 2 

2 

1 




2 

2 3 

3 


-«3 

+6 


(40) 


V. Conclusions 

In this first report we have outlined a simple, analytic 
method for obtaining the temperature distribution in a typical 
arrangement of antenna structural back*up member. The 
results were verified by the numerical shooting method. Good 
agreement between the two methods was obtained. Further 
work wfll deal with setting up the matrices for more complex 
structures, with radiation exchange, and with comparison 
between simulations and actual field measurements. 
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TaMal. Example (rflinipenrturadtatrilMitkin Table2. Examptoollsmparatimdteiributlan 

obtained by shooting mattiod (B) using tbe iinaartzation method (C) 


JT 

T 

X 


T 

0.( M ) 

0.56( H 03 

0.00 


0.560+03 

0 10 

0 559403 

0.10 


0.559+03 

0.2<» 

0.558403 

0.20 


0 558+03 

0 30 

0.557+< j 3 

0.30 


0.557+ fl 3 

0.40 

0.557403 

0.40 


0.556+03 

050 

0.556403 

0 50 


0.556+03 

0.60 

0.556403 

0.60 


0.555+03 

0.70 

0.555403 

070 


0.555+03 

0.80 

0.555+03 

0.80 


555^03 

090 

0 555+03 

090 


0.555+03 

i 00 

0.555403 

1.00 


0.554+03 

MO 

0.555403 

1 10 


0.554+03 

1.20 

O . SS 5+03 

1 20 


0.554+03 

1 30 

0 555+03 

1 30 


0 555403 

1 40 

0.55 S +(»3 

1 40 


0.555+03 

1.50 

0 555+03 

1 50 


0.555+03 

1.60 

0.555+<»3 

1 60 


0 555+03 

1 70 

0.556+03 

1.70 


0.556+03 

1 80 

0.556+03 

1 80 


0.556+03 

1 90 

0.557+03 

1 90 


0.557^)3 

2.00 

0.55 R +03 

200 


0.558+03 

2 10 

0.559+03 

2 10 


0 558+03 

2 20 

0.560+03 

2.20 


0.560+03 

2.30 

0.561+03 

2 30 


0.561+03 

2.40 

0.562+03 

2 40 


0.562+03 

2 50 

0.564+03 

2.50 


0.5649 03 

260 

0.566+03 

2 60 


0 566+03 

2 70 

0 566+03 

2 70 


0.569+03 

2 80 

0.572+03 

2.80 


0.572+03 

2 

0.576+03 

2 90 


0.576+03 

3-00 

0.580+03 

3.00 


0.580+03 
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Fig. 1. Physical system for a simple bar 


LENGTH^ m 


0.15 0.30 0.46 0.61 0.76 0.91 



Fig. 2. Temperature distribution in a simple bar using 
linearization method 



Fig. 3. Three-bar linkage 
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